defaultEncoding <- "UTF8"
library(dplyr)
library(sf)
library(terra)
library(raster)
library(DT)
library(ggplot2)
library(plotly)
library(leaflet)
casos <-
  st_read(
    "/vsicurl/https://marcelocubero.github.io/capas_proyecto/casos.geojson",
    quiet = TRUE
  )
casos <-
  casos %>%
  st_transform(4326)

casos_ <-
    read.csv( file= "https://marcelocubero.github.io/capas_proyecto/casos_.csv",
              sep = ",", encoding = "UTF-8") 


fallecidos <- rast("/vsicurl/https://marcelocubero.github.io/capas_proyecto/fallecidos.tif")
fallecidos_rl <- raster::raster(fallecidos)
bins <- c(10, 100, 500, 1000, 4000, 7600)
pal <- colorBin("Blues", domain = casos$activos, bins = bins)

leaflet() %>%
  addTiles(group="OSM") %>% 
  addProviderTiles(providers$Esri.NatGeoWorldMap , group="NatGeo") %>% 
  addProviderTiles(providers$CartoDB.DarkMatter, group="CartoDB-Black") %>%
  addRasterImage(
    fallecidos_rl, 
    opacity = 0.7,
    group= "Fallecidos",
  ) %>%
  addPolygons(
    data = casos,
    color = "black",
    fillColor = ~pal(activos),
    fillOpacity = 0.7,
    weight = 1,
    opacity = 1,
    stroke = TRUE,
    group= "Casos Activos"
  ) %>%
  addLegend(pal = pal, values = casos$activos, opacity = 0.7, title = "Casos Activos") %>%
  addLayersControl(
    baseGroups = c("OSM", "NatGeo", "CartoDB-Black"),
    overlayGroups = c("Fallecidos", "Casos Activos"),
    options = layersControlOptions(collapsed = FALSE)
  ) 
casos_ %>%
  select(Cantón= canton,Provincia= provincia, Casos_Activos= activos, Casos_Positivos= positivos) %>%
datatable( casos_, options = list(
  language = list(url = '//cdn.datatables.net/plug-ins/1.10.11/i18n/Spanish.json'),
  pageLength = 5
))
ggplot(casos_, aes(x= canton, y= positivos)) +
  geom_col( width = 0.5, fill= "#008080") +
  ggtitle("Casos Positivos por cantón") +
  xlab("Cantón") +
  ylab("Cantidad de casos positivos") +
   theme(
    plot.title = element_text(size = 16, face = "bold"),
    text = element_text(size=8),
    axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1)
  ) 

ggplot(casos_, aes(x= canton, y= activos)) +
  geom_col( width = 0.5, fill= "#00FF7F") +
  ggtitle("Casos Activos por cantón") +
  xlab("Cantón") +
  ylab("Cantidad de casos activos") +
   theme(
    plot.title = element_text(size = 16, face = "bold"),
    text = element_text(size=8),
    axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1)
  ) 

bins <- c(10, 100, 500, 1000, 4000, 7600)
pal <- colorBin("Blues", domain = casos$activos, bins = bins)


leaflet() %>%
  addTiles() %>%
  addPolygons(
    data = casos,
    color = "black",
    fillColor = ~pal(activos),
    fillOpacity = 0.7,
    weight = 1,
    opacity = 1,
    stroke = TRUE,
    group= "Casos"
  ) %>%
  addLegend(pal = pal, values = casos$activos, opacity = 0.7, title = "Casos Activos")
at <- seq(1:800)
pal1 <- colorBin('Accent', domain = at ,na.color = "transparent")

 
leaflet() %>%
  addTiles() %>% 
  addRasterImage(
    fallecidos_rl, 
    opacity = 1,
    group= "Fallecidos",
    colors = pal1
  ) %>%
  addLegend("bottomleft", pal = pal1, values = at, title = "Fallecidos")
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137
## +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null
## +wktext +no_defs +type=crs
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum World Geodetic System 1984 in Proj4 definition
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137
## +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null
## +wktext +no_defs +type=crs
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum World Geodetic System 1984 in Proj4 definition
fallecidos_rl <- raster::raster(fallecidos)
bins <- c(10, 100, 500, 1000, 4000, 7600)
pal <- colorBin("Blues", domain = casos$activos, bins = bins)
at <- seq(1:800)
pal1 <- colorBin('Accent', domain = at ,na.color = "transparent")


leaflet() %>%
  addTiles(group="OSM") %>% 
  addProviderTiles(providers$Esri.NatGeoWorldMap , group="NatGeo") %>% 
  addProviderTiles(providers$CartoDB.DarkMatter, group="CartoDB-Black") %>%
  addRasterImage(
    fallecidos_rl, 
    opacity = 1,
    group= "Fallecidos",
    colors = pal1
  ) %>%
  addLegend("bottomleft", pal = pal1, values = at, title = "Fallecidos") %>%
  addPolygons(
    data = casos,
    color = "black",
    fillColor = ~pal(activos),
    fillOpacity = 0.7,
    weight = 1,
    opacity = 1,
    stroke = TRUE,
    group= "Casos Activos"
  ) %>%
  addLegend(pal = pal, values = casos$activos, opacity = 0.7, title = "Casos Activos") %>%
  addLayersControl(
    baseGroups = c("OSM", "NatGeo", "CartoDB-Black"),
    overlayGroups = c("Fallecidos", "Casos Activos"),
    options = layersControlOptions(collapsed = FALSE)
  ) %>%
  addScaleBar("bottomright")
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137
## +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null
## +wktext +no_defs +type=crs
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum World Geodetic System 1984 in Proj4 definition
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137
## +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null
## +wktext +no_defs +type=crs
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum World Geodetic System 1984 in Proj4 definition